Genomic determinants of antibody response to a typhoid vaccine in Indian recipients.

Typhoid is endemic in India and has high global incidence. There were large outbreaks of typhoid in India between 1990 and 2018. Available typhoid vaccines induce variable levels of protective antibodies among recipients; thus, there is variability in response to the vaccine. Interindividual genomic differences is hypothesized to be a determinant of the variability in response. We studied the antibody response of ~1000 recipients of the Vi-polysaccharide typhoid vaccine from Kolkata, India, who showed considerable variability of antibody response, i.e., anti-Vi-polysaccharide antibody level 28 days postvaccination relative to prevaccination. For each vaccinee, wholegenome genotyping was performed using the Infinium Global Screening Array (Illumina). We identified 39 SNPs that mapped to 13 chromosomal regions to be associated with antibody response to the vaccine; these included SNPs on genes LRRC28 (15q26.3), RGS7 (1q43), PTPRD (9p23), CERKL (2q31.3), DGKB (7p21.2), and TCF4 (18q21.2). Many of these loci are known to be associated with various blood cell traits, autoimmune traits and responses to other vaccines; these genes are involved in immune related functions, including TLR response, JAK-STAT signalling, phagocytosis and immune homeostasis.

Single-cell transcriptomic analysis of gingivo-buccal oral cancer reveals two dominant cellular programs.

Oral squamous cell carcinoma of the gingivo-buccal region (OSCC-GB) is the most common cancer among men in India, and is associated with poor prognosis and frequent recurrence. Cellular heterogeneity in OSCC-GB was investigated by single-cell RNA sequencing of tumors derived from the oral cavity of 12 OSCC-GB patients, 3 of whom had concomitant presence of a precancerous lesion (oral submucous fibrosis [OSMF]). Unique malignant cell types, features, and phenotypic shifts in the stromal cell population were identified in oral tumors with associated submucous fibrosis. Expression levels of FOS, ATP1A, and DUSP1 provided robust discrimination between tumors with or without the concomitant presence of OSMF. Malignant cell populations shared between tumors with and without OSMF were enriched with the expression of partial epithelial-mesenchymal transition (pEMT) or fetal cell type signatures indicative of two dominant cellular programs in OSCC-GB-pEMT and fetal cellular reprogramming. Malignant cells exhibiting fetal cellular and pEMT programs were enriched with the expression of immune-related pathway genes known to be involved in antitumor immune response. In the tumor microenvironment, higher infiltration of immune cells than the stromal cells was observed. The T cell population was large in tumors and diverse subtypes of T cells with varying levels of infiltration were found. We also detected double-negative PLCG2+ T cells and cells with intermediate M1-M2 macrophage polarization. Our findings shed light on unique aspects of cellular heterogeneity and cell states in OSCC-GB.

Genomic associations with antibody response to an oral cholera vaccine.

Oral cholera vaccine is one of the key interventions used in our fight to end the longest pandemic of our time, cholera. The immune response conferred by the currently available cholera vaccines, as measured by serum antibody levels, is variable amongst its recipients. We undertook a genome wide association study (GWAS) on antibody response to the cholera vaccine; globally, the first GWAS on cholera vaccine response. We identified three clusters of bi-allelic SNPs, in high within-cluster linkage disequilibrium that were moderately (p < 5 × 10-6) associated with antibody response to the cholera vaccine and mapped to chromosomal regions 4p14, 4p16.1 and 6q23.3. Intronic SNPs of TBC1D1 comprised the cluster on 4p14, intronic SNPs of TBC1D14 comprised that on 4p16.1 and SNPs upstream of TNFAIP3 formed the cluster on 6q23.3. SNPs within and around these clusters have been implicated in immune cell function and immunological aspects of autoimmune or infectious diseases (e.g., diseases caused by Helicobacter pylori and malarial parasite). 6q23.3 is a prominent region harbouring many loci associated with immune related diseases, including multiple sclerosis, rheumatoid arthritis and systemic lupus erythematosus, as well as IL2 and INFα response to a smallpox vaccine. The gene clusters identified in this study play roles in vesicle-mediated pathway, autophagy and NF-κB signaling. No significant effect of O blood group on antibody response to the cholera vaccine was observed in this study.

Identification of HPV16 positive cervical cancer subsets characterized by divergent immune and oncogenic phenotypes with potential implications for immunotherapy.

BACKGROUND: Cervical cancers (CaCx), like many other cancer types, portray high molecular heterogeneity that affects response to therapy, including immunotherapy. In India and other developing countries, CaCx mortality rates are very high because women report to the clinics with advanced cancers in absence of organized screening programs. This calls for implementation of newer therapeutic regimens for CaCx, like immunotherapy, which is again not used commonly in such countries. OBJECTIVE: Therefore, we focused on dissecting tumour immune heterogeneity, if any, identify immune gene-based biomarkers of heterogeneity and subsets of such cancers with the potential for immunotherapy. We also attempted to characterize the cancer-associated phenotypes of such subsets, including viral load, to decipher the relationship of tumour immunogenicity with oncogenicity. METHODS: Employing RNA-seq analysis of 44 HPV16 positive CaCx patients, immune subtypes were identified by unsupervised hierarchical clustering of global immune-gene expression profiles. Proportions of tumor infiltrating immune cells in the tumor milieu were estimated, employing Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), using gene expression data from RNA-seq. The oncogenic phenotypes of the immune subtypes of CaCx were deciphered through differential gene expression (DEGs) and pathway enrichment analysis. Viral load was estimated through TaqMan-based qRT-PCR analysis. RESULTS: Analysis revealed the presence of two immune subtypes of CaCx, A (26/44; 59.09%) and B (18/44; 40.90%). Compared to Subtype-A, Subtype-B portrayed overexpression of immune genes and high infiltration of immune cells, specifically CD8+ T cells (p < 0.0001). Besides, a significant correlation between PD-1 and PD-L1 co-expression among Subtype-B, as opposed to Subtype-A, confirmed the interactive roles of these immune checkpoint molecules in Subtype B. Stepwise discriminant analysis pin-pointed ten immune-genes that could classify 100% of the patients significantly (p < 0.0001) into the two immune subtypes and serve as potential biomarkers of CaCx immunity. Differential gene expression analysis between the subtypes unveiled that Subtype-B was more biologically aggressive than Subtype-A, reflecting loss of structural integrity and promotion of cancer progression. The viral load was significantly lower in Subtype-B (average viral load = 10.74/100 ng of genomic DNA) compared to Subtype-A (average viral load = 14.29/100 ng of genomic DNA). Thus viral load and the ten-gene panel underscore their association with immunogenicity and oncogenicity. CONCLUSION: Our study provides strong evidence that only a subset, about 41% of HPV16 positive CaCx patients in India, portray immune enrichment of the tumor milieu coupled with aggressive phenotypes. Such subtypes are therefore likely to benefit through checkpoint molecule-based or tumor infiltrating lymphocyte-based immunotherapy, which could be a leap forward in tackling aggressive forms of such CaCx in India and other developing countries.

Overexpression of CD73 is associated with recurrence and poor prognosis of gingivobuccal oral cancer as revealed by transcriptome and deep immune profiling of paired tumor and margin tissues.

BACKGROUND: For various cancers, differences in response to treatment and subsequent survival period have been reported to be associated with variation in immune contextures. AIM: We sought to identify whether such association exists in respect of gingivobuccal oral cancer. MATERIALS AND METHODS: We performed deep immune profiling of tumor and margin tissues collected from 46 treatment naïve, Human Papillomavirus (HPV) negative, patients. Each patient was followed for 24 months and prognosis (recurrence/death) noted. Key findings were validated by comparing with TCGA-HNSC cohort data. RESULTS: About 28% of patients showed poor post-treatment prognosis. These patients exhibited a high probability of recurrence even within 1 year and death within 2 years. There was restricted immune cell infiltration in tumor, but not in margin, among these patients. Reduced expression of eight immune-related genes (IRGs) (NT5E, THRA, RBP1, TLR4, ITGA6, BMPR1B, ITGAV, SSTR1) in tumor strongly predicted better quality of prognosis, both in our patient cohort and in TCGA-HNSC cohort. Tumors of patients with better prognosis were associated with (a) lower CD73+ cells with concomitant lower expression level of NT5E/CD73, (b) higher proportions of CD4+ and CD8+ T cells, B cells, NK cells, M1 macrophages, (c) higher %Granzyme+ cells, (d) higher TCR and BCR repertoire diversities. CD73 expression in tumor was associated with low CD8+ and CD4+ T cells, low immune repertoire diversity, and advanced cancer stage. DISCUSSION AND CONCLUSION: High infiltration of anti-tumor immune cells in both tumors and margins results in good prognosis, while in patients with minimal infiltration in tumors in spite of high infiltration in margins results in poor prognosis. Targeted CD73 immune-checkpoint inhibition may improve clinical outcome.

Genetic variants associated with spontaneous preterm birth in women from India: a prospective cohort study.

Background: Despite having the highest number of preterm births globally, no genomic study on preterm birth was previously published from India or other South-Asian countries. Methods: We conducted a genome-wide association (GWA) study of spontaneous preterm birth (sPTB) on 6211 women from India. We used a novel resampling procedure to identify the associated single nucleotide polymorphisms (SNPs) followed by haplotype association analysis and imputation. Findings: We found that 512 maternal SNPs were associated with sPTB (p < 2.51e-3), of which minor allele at 19 SNPs (after Bonferroni correction) had increased genotype relative risk. Haplotypes containing six of the 19 SNPs (rs13011430, rs8179838, rs2327290, rs4798499, rs7629800, and rs13180906) were associated with sPTB (p < 9.9e-4; Bonferroni adjusted p-value <0.05). After imputation in regions around the 19 SNPs, 15 imputed SNPs were found to be associated with sPTB (Bonferroni adjusted p-value <0.05). One of these imputed SNPs, rs35760881, and three other SNPs (rs17307697, rs4308815, and rs10983507) were also reported to be associated with sPTB in women belonging to European ancestry. Moreover, we found that GG genotype at rs1152954, one of the associated SNPs, enhanced risk of sPTB and reduced telomere length. Interpretation: This is the first study from South Asia on the genome-wide identification of maternal SNPs associated with sPTB. These SNPs are known to alter the expression of genes associated with major pathways in sPTB viz. inflammation, apoptosis, cervical ripening, telomere maintenance, selenocysteine biosynthesis, myometrial contraction, and innate immunity. From a public health perspective, the trans-ethnic association of four SNPs identified in our study may help to stratify women with risk of sPTB in most populations. Funding: Department of Biotechnology (India), Grand Challenges India - All Children Thriving Program and Biotechnology Industry Research Assistance Council (BIRAC).

Genome-wide temporal landscaping of DNA methylation in pregnant women delivering at term: a GARBH-InI study.

Background: We performed an epigenome-wide longitudinal DNA methylation study on an Indian cohort of pregnant women, GARBH-Ini, at three time points during pregnancy and at delivery. Aim & objective: Our aim was to identify temporal DNA methylation changes in maternal peripheral blood during the period of gestation and assess their impact on biological pathways critical for term delivery. Results: Significantly differentially methylated CpGs were identified by linear mixed model analysis (Bonferroni p < 0.01) and classified into two distinct temporal methylation trends: increasing and decreasing during gestation. Genes with upward methylation trend were enriched for T-cell activity, while those with a downward trend were enriched for solute transport and cell structure organization functions. Conclusion: Consistent trends of DNA methylation in maternal peripheral blood point to the sentinel function of T cells in the maintenance of pregnancy, and the importance of coordinated cellular remodeling to facilitate term delivery.

DNA methylation is the addition of a methyl group to the molecular structure of DNA, which then alters the gene expression. The goal of the study was to find out how DNA methylation patterns change over time during pregnancy and how these changes are related to the biological processes that are important for the delivery of a healthy baby at full term. Using statistical modeling, we identified specific patterns of DNA methylation changes during pregnancy and classified them into two groups based on the direction of the changes. The genes associated with increasing methylation levels were related to the activities of T cells, which are important for the immune system. The genes associated with decreasing methylation levels were related to processes like transporting substances and organizing cell structures. In conclusion, our findings suggest that T cells play an important role in maintaining a healthy pregnancy, and the study highlights the importance of coordinated changes in cells to support a successful delivery of a baby at term.

Two-step approach based multi-objective groundwater remediation using enhanced random vector functional link integrated with evolutionary marine predator algorithm.

We here propose a two-step approach-based simulation-optimization model for multi-objective groundwater remediation using enhanced random vector functional link (ERVFL) and evolutionary marine predator algorithm (EMPA). In this study, groundwater flow and solute transport models are developed using MODFLOW and MT3DMS. The ERVFL network is used to approximate the flow and transport models, enhancing the computational performance. This study also improves the robustness of the ERVFL network using a kernel density estimator (KDE) based weighted least square approach. We further develop the EMPA by modifying the marine predator algorithm (MPA) using elite opposition-based learning, biological evolution operators, and elimination mechanisms. In the multi-objective version of EMPA, the non-dominated/Pareto-optimal solutions are stored in an external repository using an archive controller and adaptive grid mechanism to promote better convergence and diversity of the Pareto front. The proposed methodologies are applied for multi-objective groundwater remediation of a hypothetical unconfined aquifer based on the two-step method. The first step directly integrates flow and transport models with EMPA and finds the optimal locations of pumping wells by minimizing the percent of contaminant mass remaining in the aquifer. In the second step, the ERVL-based proxy model is integrated with EMPA and used for multi-objective optimization while explicitly using the pumping well locations obtained in the first step. The multi-objective optimization generates a Pareto-optimal solution representing the relationship between the rate of pumping and the amount of contaminant mass in the aquifer. Further analyses show a significant advantage of the two-step approach over a traditional method for multi-objective groundwater remediation.

Twitterati on COVID-19 pandemic-environment linkage: Insights from mining one year of tweets.

The COVID-19 pandemic seems to have had positive (although short-lived, e.g., reduction in pollution due to lockdown) as well as negative (e.g., increasing plastic pollution due to use of disposable masks, etc.) impacts on the environment. The pandemic-environment linkage also includes circumstances when regions experienced extreme weather events, such as floods and cyclones, and disaster management became challenging. This study aims to examine the trends in public discourses on Twitter on these interactions between the pandemic and environment. The present study follows the most recent literature on understanding public perceptions - which acknowledges Twitter to be an abundant source of information on public discussions on any global issue, including the pandemic. A Python-based code is developed to extract Twitter data spanning over a year, and analyze the presence of covid-environment related keywords and other attributes. It is found that the Twitterati aggressively viewed the impacts (such as economic slowdown and high mortality) of the pandemic as miniatures of the results of future climate change. The community was also highly concerned about the varying air and plastic pollution levels with the change in lockdown and covid prevention policies. Extreme weather events were a high-frequency topic when they impacted countries such as India, the USA, Australia, the Philippines and Vietnam. This study makes a novel attempt to provide an overview of public discourses on the pandemic-environment linkage and; can be a crucial addition to the literature on assessing public perception of environmental threats through Twitter data mining.

Multiple migrations from East Asia led to linguistic transformation in NorthEast India and mainland Southeast Asia.

NorthEast India, with its unique geographic location in the midst of the Himalayas and Bay of Bengal, has served as a passage for the movement of modern humans across the Indian subcontinent and East/Southeast Asia. In this study we look into the population genetics of a unique population called the Khasi, speaking a language (also known as the Khasi language) belonging to the Austroasiatic language family and residing amidst the Tibeto-Burman speakers as an isolated population. The Khasi language belongs to one of the three major broad classifications or phyla of the Austroasiatic language and the speakers of the three sub-groups are separated from each other by large geographical distances. The Khasi speakers are separated from their nearest Austroasiatic language-speaking sub-groups: the "Mundari" sub-family from East and peninsular India and the "Mon-Khmers" in Mainland Southeast Asia. We found the Khasi population to be genetically distinct from other Austroasiatic speakers, i.e. Mundaris and Mon-Khmers, but relatively similar to the geographically proximal Tibeto Burmans. The possible reasons for this genetic-linguistic discordance lie in the admixture history of different migration events that originated from East Asia and proceeded possibly towards Southeast Asia. We found at least two distinct migration events from East Asia. While the ancestors of today's Tibeto-Burman speakers were affected by both, the ancestors of Khasis were insulated from the second migration event. Correlating the linguistic similarity of Tibeto-Burman and Sino-Tibetan languages of today's East Asians, we infer that the second wave of migration resulted in a linguistic transition while the Khasis could preserve their linguistic identity.

Integrative analysis of genomic and transcriptomic data of normal, tumour, and co-occurring leukoplakia tissue triads drawn from patients with gingivobuccal oral cancer identifies signatures of tumour initiation and progression.

A thickened, white patch - leukoplakia - in the oral cavity is usually benign, but sometimes (in ~9% of individuals) it progresses to malignant tumour. Because the genomic basis of this progression is poorly understood, we undertook this study and collected samples of four tissues - leukoplakia, tumour, adjacent normal, and blood - from each of 28 patients suffering from gingivobuccal oral cancer. We performed multiomics analysis of the 112 collected tissues (four tissues per patient from 28 patients) and integrated information on progressive changes in the mutational and transcriptional profiles of each patient to create this genomic narrative. Additionally, we generated and analysed whole-exome sequence data from leukoplakia tissues collected from 11 individuals not suffering from oral cancer. Nonsynonymous somatic mutations in the CASP8 gene were identified as the likely events to initiate malignant transformation, since these were frequently shared between tumour and co-occurring leukoplakia. CASP8 alterations were also shown to enhance expressions of genes that favour lateral spread of mutant cells. During malignant transformation, additional pathogenic mutations are acquired in key genes (TP53, NOTCH1, HRAS) (41% of patients); chromosomal-instability (arm-level deletions of 19p and q, focal-deletion of DNA-repair pathway genes and NOTCH1, amplification of EGFR) (77%), and increased APOBEC-activity (23%) are also observed. These additional alterations were present singly (18% of patients) or in combination (68%). Some of these alterations likely impact immune-dynamics of the evolving transformed tissue; progression to malignancy is associated with immune suppression through infiltration of regulatory T-cells (56%), depletion of cytotoxic T-cells (68%), and antigen-presenting dendritic cells (72%), with a concomitant increase in inflammation (92%). Patients can be grouped into three clusters by the estimated time to development of cancer from precancer by acquiring additional mutations (range: 4-10 years). Our findings provide deep molecular insights into the evolutionary processes and trajectories of oral cancer initiation and progression. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

A novel two-step approach for optimal groundwater remediation by coupling extreme learning machine with evolutionary hunting strategy based metaheuristics.

We propose a simulation-optimization (SO) model based on a novel two-step strategy for the optimal design of groundwater remediation systems. The SO models are developed by coupling simulation models directly or through the extreme learning machine (ELM) with evolutionary hunting strategy based metaheuristics (EHSMs). In the first step, EHSMs with a combinatorial optimization technique are used to obtain optimal pumping locations by minimizing the percentage of contaminant mass that remained in the aquifer while keeping the pumping strategy as constant. In the second step, the optimal pumping locations are directly used as input, and a composite function is employed to minimize the sum of the water extraction rates and the percentage of extracted contaminant mass by constraining hydraulic heads and contaminant concentrations. The performance of the two-step strategy is found to be slightly better and computationally more efficient than the alternate approach. Moreover, various statistical measures suggest the superiority of EHSMs over other metaheuristics for groundwater remediation.

Genes and pathways monotonically dysregulated during progression from normal through leukoplakia to gingivo-buccal oral cancer.

Oral squamous cell carcinoma of the gingivo-buccal region (OSCC-GB) accounts for the highest cancer morbidity and mortality among men in India. It has been observed that about one-third of individuals with oral leukoplakia, a dysplastic precancerous lesion in the oral cavity, progress to oral cancer. We aimed to identify systematic transcriptomic changes as a normal tissue in the oral cavity progresses to frank OSCC-GB. Seventy-two OSCC-GB patients, from multiple hospitals, were recruited, and transcriptome analysis of tumor and adjacent normal tissue (of all patients) and adjacent leukoplakia tissue (of a subset of 25 unselected patients with concomitant leukoplakia) was performed. We have identified many differences in the transcriptomic profiles between OSCC-GB and squamous cell carcinoma of the head and neck regions. Compared to the normal/precancerous tissue, significant enrichment of ECM-receptor interaction, PI3K-Akt signaling, cytokine-cytokine receptor interaction, focal adhesion, and cell cycle pathways were observed in OSCC-GB. Using gene set enrichment analysis, we identified a profound role of interferon receptor signaling in tumor growth by activating immune evasion mechanisms. The role of tumor-infiltrating immune cells further supported the growth and immunosuppressive mechanism of tumor tissues. Some immune evasion genes-CD274, CD80, and IDO1-were found to be activated even in the precancerous tissue. Taken together, our findings provide a clear insight into the sequential genetic dysregulation associated with progression to oral cancer. This insight provides a window to the development of predictive biomarkers and therapeutic targets for gingivo-buccal oral cancer.

dbGENVOC: database of GENomic Variants of Oral Cancer, with special reference to India.

Oral cancer is highly prevalent in India and is the most frequent cancer type among Indian males. It is also very common in southeast Asia. India has participated in the International Cancer Genome Consortium (ICGC) and some national initiatives to generate large-scale genomic data on oral cancer patients and analyze to identify associations and systematically catalog the associated variants. We have now created an open, web-accessible database of these variants found significantly associated with Indian oral cancer patients, with a user-friendly interface to enable easy mining. We have value added to this database by including relevant data collated from various sources on other global populations, thereby providing opportunities of comparative geographical and/or ethnic analyses. Currently, no other database of similar nature is available on oral cancer. We have developed Database of GENomic Variants of Oral Cancer, a browsable online database framework for storage, retrieval and analysis of large-scale data on genomic variants and make it freely accessible to the scientific community. Presently, the web-accessible database allows potential users to mine data on ∼24 million clinically relevant somatic and germline variants derived from exomes (n = 100) and whole genomes (n = 5) of Indian oral cancer patients; all generated by us. Variant data from The Cancer Genome Atlas and data manually curated from peer-reviewed publications were also incorporated into the database for comparative analyses. It allows users to query the database by a single gene, multiple genes, multiple variant sites, genomic region, patient ID and pathway identities. Database URL: http://research.nibmg.ac.in/dbcares/dbgenvoc/.

SARS-CoV-2 mutation 614G creates an elastase cleavage site enhancing its spread in high AAT-deficient regions.

SARS-CoV-2 was first reported from China. Within three months, it evolved to 10 additional subtypes. Two evolved subtypes (A2 and A2a) carry a non-synonymous Spike protein mutation (D614G). We conducted phylodynamic analysis of over 70,000 SARS-CoV-2 coronaviruses worldwide, sequenced until July2020, and found that the mutant subtype (614G) outcompeted the pre-existing type (614D), significantly faster in Europe and North-America than in East Asia. Bioinformatically and computationally, we identified a novel neutrophil elastase (ELANE) cleavage site introduced in the G-mutant, near the S1-S2 junction of the Spike protein. We hypothesised that elevation of neutrophil elastase level at the site of infection will enhance the activation of Spike protein thus facilitating host cell entry for 614G, but not the 614D, subtype. The level of neutrophil elastase in the lung is modulated by its inhibitor α1-antitrypsin (AAT). AAT prevents lung tissue damage by elastase. However, many individuals exhibit genotype-dependent deficiency of AAT. AAT deficiency eases host-cell entry of the 614G virus, by retarding inhibition of neutrophil elastase and consequently enhancing activation of the Spike protein. AAT deficiency is highly prevalent in European and North-American populations, but much less so in East Asia. Therefore, the 614G subtype is able to infect and spread more easily in populations of the former regions than in the latter region. Our analyses provide a molecular biological and evolutionary model for the higher observed virulence of the 614G subtype, in terms of causing higher morbidity in the host (higher infectivity and higher viral load), than the non-mutant 614D subtype.

Application of Random Forest and data integration identifies three dysregulated genes and enrichment of Central Carbon Metabolism pathway in Oral Cancer.

BACKGROUND: Studies of epigenomic alterations associated with diseases primarily focus on methylation profiles of promoter regions of genes, but not of other genomic regions. In our past work (Das et al. 2019) on patients suffering from gingivo-buccal oral cancer - the most prevalent form of cancer among males in India - we have also focused on promoter methylation changes and resultant impact on transcription profiles. Here, we have investigated alterations in non-promoter (gene-body) methylation profiles and have carried out an integrative analysis of gene-body methylation and transcriptomic data of oral cancer patients. METHODS: Tumor and adjacent normal tissue samples were collected from 40 patients. Data on methylation in the non-promoter (gene-body) regions of genes and transcriptome profiles were generated and analyzed. Because of high dimensionality and highly correlated nature of these data, we have used Random Forest (RF) and other data-analytical methods. RESULTS: Integrative analysis of non-promoter methylation and transcriptome data revealed significant methylation-driven alterations in some genes that also significantly impact on their transcription levels. These changes result in enrichment of the Central Carbon Metabolism (CCM) pathway, primarily by dysregulation of (a) NTRK3, which plays a dual role as an oncogene and a tumor suppressor; (b) SLC7A5 (LAT1) which is a transporter dedicated to essential amino acids, and is overexpressed in cancer cells to meet the increased demand for nutrients that include glucose and essential amino acids; and, (c) EGFR which has been earlier implicated in progression, recurrence, and stemness of oral cancer, but we provide evidence of epigenetic impact on overexpression of this gene for the first time. CONCLUSIONS: In rapidly dividing cancer cells, metabolic reprogramming from normal cells takes place to enable enhanced proliferation. Here, we have identified that among oral cancer patients, genes in the CCM pathway - that plays a fundamental role in metabolic reprogramming - are significantly dysregulated because of perturbation of methylation in non-promoter regions of the genome. This result compliments our previous result that perturbation of promoter methylation results in significant changes in key genes that regulate the feedback process of DNA methylation for the maintenance of normal cell division.

Dysbiotic Lesional Microbiome With Filaggrin Missense Variants Associate With Atopic Dermatitis in India.

Background: Atopic Dermatitis (AD) has been associated with the loss of function (LoF) mutations in Filaggrin (FLG) gene and increase in relative abundance of specific microbes in the lesional skin, predominantly in Caucasians. Our study aims to determine, in Indian AD patients, (a) the prevalence of FLG LoF and missense mutations, and (b) the nature and extent of dysbiosis and altered microbial pathways with and without mutations in FLG. AD patients (n = 34) and healthy controls (n = 54) were recruited from India in this study and shotgun sequencing was carried out in a subset of samples with adequate microbiome DNA concentration. Host DNA from the same subset of samples was subjected to FLG coding region sequencing and host-microbiome association was estimated. Results: The prevalence of FLG LoFs that are associated with AD globally were significantly lesser in our cases and controls (8.6%, 0%) than those reported in Europeans (27%, 2.6%). Staphylococcus aureus was present only on AD skin [abundance in Pediatric AD: 32.86%; Adult AD: 22.17%], but not on healthy skin on which Staphylococcus hominis (Adult controls: 16.43%, Adult AD: 0.20%; p = 0.002), Cutibacterium acnes (Adult controls:10.84%, Adult AD: 0.90%; p = 0.02), and Malassezia globosa (Adult controls: 8.89%, Adult AD: 0.005%; p = 0.001) were significantly more abundant. Microbial pathways mostly associated with skin barrier permeability, ammonia production and inflammation (Arginine and Proline metabolism, Histidine Metabolism and Staphylococcus aureus infection) were significantly enriched on AD skin metagenome. These pathways are also reported to impair antimicrobial peptide activity. Among AD patients with missense single nucleotide polymorphisms harboring "potentially damaging" alleles in FLG gene, damaging allele dosage was significantly (p < 0.02) positively correlated with relative abundance of phylum_Proteobacteria up to order_Pseudomonadales and negatively correlated with phylum_Firmicutes up to species_Staphylococcus aureus. Conclusion: Our study has provided evidence that host DNA profile is significantly associated with microbiome composition in the development of AD. Species and strain level analysis showed that the microbial pathways enriched in AD cases were mostly found in MRSA strains. These evidences can be harnessed to control AD by modulating the microbiome using a personalized strategy. Our findings on the association of FLG genotypes with the microbiome dysbiosis may pave the way for a personalized strategy to provide a more effective control of AD.

Density of CD3+ and CD8+ cells in gingivo-buccal oral squamous cell carcinoma is associated with lymph node metastases and survival.

The tumor immune microenvironment is emerging as a critical player in predicting cancer prognosis and response to therapies. However, the prognostic value of tumor-infiltrating immune cells in Gingivo-Buccal Oral Squamous Cell Carcinoma (GBOSCC) and their association with tumor size or lymph node metastases status require further elucidation. To study the relationship of tumor-infiltrating immune cells with tumor size (T stage) and lymph node metastases (N stages), we analyzed the density of tumor-infiltrating immune cells in archived, whole tumor resections from 94 patients. We characterized these sections by immune-histochemistry using 12 markers and enumerated tumor-infiltrating immune cells at the invasive margins (IM) and centers of tumors (CT). We observed that a higher density of CD3+ cells in the IM and CT was associated with smaller tumor size (T1-T2 stage). Fewer CD3+ cells was associated with larger tumor size (T3-T4 stage). High infiltration of CD3+and CD8+ cells in IM and CT as well as high CD4+ cell infiltrates in the IM was significantly associated with the absence of lymph node metastases. High infiltrates of CD3+ and CD8+ cells in CT was associated with significantly improved survival. Our results illustrate that the densities and spatial distribution of CD3+ and CD8+ cell infiltrates in primary GBOSCC tumors is predictive of disease progression and survival. Based on our findings, we recommend incorporating immune cell quantification in the TNM classification and routine histopathology reporting of GBOSCC. Immune cell quantification in CT and IM may help predict the efficacy of future therapies.